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Abstract. We report the identification and construction of raising and lowering operators 
for the complete eigenfunctions of isotropic harmonic oscillators confined by dihedral angles, 
in circular cylindrical and spherical coordinates; as well as for the hydrogen atom in the same 
situation of confinement, in spherical, parabolic and prolate spheroidal coordinates. The 
actions of such operators on any eigenfunction are examined in the respective coordinates, 
illustrating the possibility of generating the complete bases of eigenfunctions in the respective 
coordinates for both physical systems. The relationships between the eigenfunctions in each 
pair of coordinates, and with the same eigenenergies are also illustrated. 
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1 Introduction 

In this section, we provide the background information on our investigations on the hydrogen 
atom in diverse situations of confinement motivated by a variety of physical and chemical situa- 
tions, spanning a period of thirty-some years as reviewed in [15], and also our current work for 
confinement by a dihedral angle [23]. In each case the Schrodinger equation is separable in the 
chosen coordinate system, and its solutions differ from the familiar ones for the free hydrogen 
atom by the condition of vanishing at the position of the confining natural boundary, determined 
by the fixed value of one of the coordinates. Such "natural" boundaries allow for separable 
and integrable solutions, with non-integer eigenvalues for the constant of motion canonically 
conjugate to the confining coordinate. The presence of such eigenvalues in the differential 
equations for the other two degrees of freedom affects correspondingly their eigensolutions, which 
nevertheless maintain their respective integer quantum numbers and constants of motion. The 
following paragraphs contain a description of our motivation to search, identify and construct 
ladder operators for both the harmonic oscillator and the hydrogen atom confined by dihedral 
angles, in the different coordinate systems sharing the broken rotational symmetry around the 
edge of the angle. 

The research reported in this article was motivated by the feedback on our contribution in 
the Symposium on Superintegrability, Exact Solvability, and Special Functions. In fact, that 
contribution with the title "Superintegrability in Confined Quantum Systems" was focused on 
the case of the hydrogen atom confined inside spheres, and prolate spheroids, and in semi- 
infinite spaces with paraboloidal, hyperboloidal and circular cone boundaries [15] as well as 

*This paper is a contribution to the Special Issue "Superintegrability, Exact Solvability, and Special Functions" . 
The full collection is available at http://www.cmis.de/journals/SIGMA/SESSF2012.html 
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dihedral angles [23]. The successive boundary condition equations determining the eigenvalues 
of the respective constants of motion canonically conjugate to the confining coordinate are 



M(-u r ,2£ + 2,2r /v) = 0, v = v r +£ + l, (1) 

S Uu (uq) = 0, v = v u + n v +m + l, (2) 

M(-vz,m + l,$/v) =0, v = u^ + n v + m + l, (3) 

Sv v (v ) = 0, v = n u + v v + m + l, (4) 

Pf(costf ) = 0, v = n r + \ + l, (5) 

sin/i^o = 0, v = n r + ng + // + ! = n„ + n„ + ^ + l = n^ + + + (6) 



involving confluent hypergeometric functions M(a,b, z), prolate spheroidal functions 5* in the 
spheroidal 1 < uq < oo and hyperboloidal — 1 ^ i>o ^ 1 coordinates, associated Legendre P™ 
functions, and trigonometric sine functions. Their respective quantum labels f r , u u , u v , A 
and fx are no longer integer. Consequently, the functions themselves are no longer polynomials 
as in the free hydrogen atom, and the principal quantum label v is no longer an integer, either. 
The energy eigenvalues E v = —e 2 /2aou 2 and their degeneracies are also changed accordingly. 

Since exact solvability requires polynomial solutions and raising and lowering operators con- 
necting them, the first five types of confinement do not belong to this category. Nevertheless, 
their solutions can be accurately evaluated for practical applications [15] and original references 
therein. 

Concerning the confinement by dihedral angles, the solutions in the other degrees of freedom 
maintain their polynomial character in spherical, spheroidal and parabolic coordinates. In order 
to establish whether they belong or not to the exactly solvable category, we must search for the 
respective raising and lowering operators. 

From the known connection of the solutions of the free hydrogen atom in parabolic coordinates 
with the solutions of the harmonic oscillator in circular coordinates as implied in equation (3), we 
have decided to search first for the corresponding operators in the isotropic harmonic oscillator 
confined by dihedral angles, in circular cylindrical and spherical coordinates. Correspondingly, 
in Section 2 we identify the eigenfunctions and energy eigenvalues in the respective coordinates. 
In Section 3 the starting point is the familiar raising and lowering operators for the free har- 
monic oscillator in cartesian coordinates, followed by the identification of their forms in circular 
cylindrical and spherical coordinates; the next step is to incorporate the restrictions associated 
with the eigenfunctions in the confining angle, and finally examine its consequences on the other 
degrees of freedom. For the latter we make full use of the results of Section 2 establishing by 
mathematical induction the explicit and general forms of the raising and lowering operators 
for the circular, spherical radial and polar angle coordinate eigenfunctions. In Section 4 the 
relationships between the sets of degenerate circular cylindrical and spherical eigenfunctions are 
explicitly illustrated. The corresponding results for the hydrogen atom confined by dihedral 
angles are presented in the successive Sections 5-7. The discussion in Section 8 revolves around 
the superintegrability, exact solvability and special functions for both confined-by-dihedral-angle 
quantum systems, and others. 

2 Eigenfunctions and eigenenergies of the isotropic 
harmonic oscillator confined by dihedral angles 

The Schrodinger equation for the harmonic oscillator is known to be separable in cartesian, 
circular cylindrical and spherical coordinates [7, 19]: 



x = p cos ip = r sin $ cos <p, y = p sin tp = r sin •& sin 99, z = r cos 



(7) 
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For the system confined by a dihedral angle defined by the meridian half-planes <p = and 
p = <po, the boundary conditions on the respective eigenfunctions are 



ip(p, <p = 0, z) = 0, tp(p, <p = <p , z) = 0, 
ip(r,0,<p = O) = 0, ^(r,^<^ = ^ ) = 0. 



(9) 



The separability of the Schrodinger equation leads to the common eigenvalue equation for the 
square of the z-component of the angular momentum 



dip'- 



(10) 



As well as the eigenvalue equations for the circular radial and longitudinal cartesian Hamiltonians 



H 2 (I d d p 2 



2m \p dp dp 

h 2 d 2 l 



l 



2 „2 



+ -muj p 



R(p) = E n ,R(p), 



(11) 



2m dz 2 2 



2 2 

H — muj z 



Z(z) = E nz Z(z), 



and the counterparts for the square of the angular momentum and spherical radial Hamiltonian, 
respectively, 



1 d 



sin 



e i + -At^ W = A ( A + !) W' 



sin^d^ ci0 sin 2 #y 
h 2 



2m 



1 d 2 d A(A + 1) 
r 2 dr dr r 2 



+ -muj 2 r 2 }R(i 1 )=E npX R(r). 



The boundary conditions of equations (8), (9) lead to the eigenfunctions and eigenvalues: 



sin pep, 



<Po 



n.. 



1,2,3,.... 



(12) 



We write first its circular radial and axial coordinate companion factors in the cylindrical case, 
with length units yfh/mtu, 



Rn P Jp) = N np ^e-P 2 ' 2 M (- np ,p + l,p 2 ) , 
Z nz {z) = N nz e~ z2 / 2 H nz {z) 



(13) 
(14) 



in terms of singularity removing factors for p — > 0, p — > oo and z — ¥ oo, and the confluent 
hypergeometric polynomial of degree 2n p and the Hermite polynomials, respectively. The latter 
correspond to the free harmonic oscillator in the cartesian axial coordinate, with parity (— ) Uz . 
The radial eigenfunction has the same form as in the free harmonic oscillator with the integer 
magnetic quantum number m replaced by the p eigenvalues of equation (12). The respective 
energy eigenvalues in units fvjj become 



E„ 



2n p + p + n z + - . 



JripUzti y • t~ • ■ 2' 

Next, we go on to the spherical radial and polar angle coordinate eigenfunctions: 



n ( r \ — at r no+ ^e~ r2 l 2 M 



3 2 

-n r ,ng + p + -,r 



(15) 



(16) 



e^ +At (cos0) = N ng+I1 



sim 



2 Fi -n e , n e + 2p + 1; p + 1; 



1 — cos I 
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Table 1. Polynomial confluent hypergeometric functions. 



n. 



M(-n p ,p + l,p 2 



1 



3p 2 



+ 



2p 2 
3P 4 



P 
M+l 



+ 



(M+l)(p+2) 



M+l ^ (/^+1)(m+2) (/i+l)(/t+2)0+3) 



Table 2. Hermite polynomials. 





Hermite(n 2 , z) 





1 


1 


2z 


2 


-2 + 4z 2 


3 


-12z + 8z 3 


4 


12 - 48z 2 + 16z 4 


5 


120z - 160z 3 + 32z 5 


6 


-120 + 720z 2 - 480z 4 + 64z 6 



+ R ne 2^1 ( -n fl) ng + 2/x + 1; fjt + 1; 



(17) 

where the identification of the angular momentum label A = ng + /x is in order. The remarks 
about the radial eigenfunction are the same as for its cylindrical counterpart of equation (13). 
Concerning the polar angle eigenfunction, several remarks are necessary: i) each of the terms 
inside the brackets corresponds to a hypergeometric function representation of the associated 
Legendre polynomials, ii) for integer values of fj, — > m they are equal to each other for ng even 
or odd, containing only even or odd powers of cos#, in) for the values of fi in equation (12), 
both contain even and odd powers of cos and do not have a well defined parity, iv) their linear 
combination in equation (17) ensures that its parity (— ) n ° is well-defined. While, the confining 
in the dihedral angle breaks the rotational symmetry around the z-axis, the /i values in each of 
the terms in equation (17) also break the z — > — z parity; restoration of such a parity requires 
the chosen superposition. The corresponding energy eigenvalues are 

3 

En r n et i = 2n r + ng + -. (18) 

For a chosen value of fj,, the states with 2n p + n z = 2n r + ng share the same energy. Notice 
that the Gaussian exponential factors in equation (13) and (14) reproduce the corresponding 
factor in equation (16); while the u powers of r and sm.0 in equations (16) and (17) reproduce 
the n power of p in equation (13). Consequently, the eigenfunctions with n p = n z = and 
n r = ng = are the same in both coordinate systems for all their common values of u, or n v 
for each confining angle (fo , equation (12). 

Tables 1-4 illustrate the polynomial factors in the lower excited radial and axial, and radial 
and polar angle polynomials equations (13), (14) and (16), (17) in the eigenfunctions for the 
respective coordinate systems. Each one of them serves as a guide in the identification of 
their respective raising and lowering operators; and the comparison of their products leads 
to the relationships between the complete sets of degenerate eigenfunctions in the alternative 
coordinates. 
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Table 3. Polynomial confluent hypergeometric functions. 



n r 


M(-n r , A + §,r 2 ) 



1 

2 

3 


1 

1- A 

A+§ 

1 2r 2 | r 4 

A+f T (A+|)(A+f) 
-i 3r 2 | 3r 4 r 6 

A+T ^ (A+|)(A+|) (A+|)(A+|)(A+|) 



Table 4. Associated Legendre polynomials with definite parity (— ) 



no 


Even or odd powers of 2-F1 (— ne, + 2/i + 1; u + 1; 1 ™ s6> ) 





1 


1 


cos # 


2 
3 
4 
5 
6 


-l+(2/^+3)cos 2 

2(M+1) 
-3 cos 6»+(2/i+5)cos 3 6 
2(^+1) 

3-6(2/x+5)cos 2 0+(2/^+5)(2/i+7)cos 4 
4(|*+l)0H-2) 
15cos0-lO(2^+7)cos- i e+(2^+7)(2/i+9)cos J 6» 
4( M +l)( M +2) 

-15+45((2^+7) cos 2 6»-15(2^+7)(2ft+9)cos 4 6l+(2/i+7)(2^+9)(2/j+ll)cos 6 


8(j^+l)(M+2)(M+3) 



3 Raising and lowering operators 

for the successive degrees of freedom 

The creation and annihilation operators for the free harmonic oscillator in each cartesian co- 
ordinate constitute our starting point and immediate bridge to their circular cylindrical and 
spherical counterparts: 

? (* ± |) + 1 {y ± |) + « (' ± |) = (? ± v,) + S (z ± |) = (f ± v„) . (19) 



C?±v,)=?(,±iL)± 



Let us consider the circular coordinate components rewritten in terms of the cartesian unit 
vectors i and j: 

d \ , tp d 

= ( i cos ip + j sin tp) ( p ± — ) ± (— i sin 03 + j cos cp) — — . (20) 
v \ op J ' pop 

The application of the operators in the i direction on the eigenfunctions of equation (12), with 
the idea of raising and lowering the quantum number n„, suggests the change of p by mp/p)Q 
in the respective trigonometric functions. In fact, for the angular factor of the radial operator 
term the result is 



nip n^irp 1 

cos — sin = - 

p p 2 



. (n v + l)np . (ritn - l)irip 
sin — h sin — 



(21) 



achieving the desired raising and lowering actions in the eigenfunctions. For the angular deriva- 
tive term, the corresponding result is similar, 

Trip d n^irp n v 7r irp n^irp 

— sin — sin — = — sin — cos — 

po dp po po <Po <Po 
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u^tt 1 
ipo 2 



. (n„, + l)7Ttf . (n<p - l)n(p 
sin — sin — 



po 



¥>0 



(22) 



The appropriate linear combinations of equations (21) and (22) allow the identification of 
the respective raising and lowering operators: 



TT(p 
COS h 



ipo Trip d 
sin 



cos 



TTip 



sm 



ip dip 
np d 
ipo dip 



n 



sm ■ 



sm 



<p7rp 
n^TTip 



sm ■ 



sm 



(n v + l)ir(p 
(n v - l)irip 



(23) 
(24) 



They have the same structure as the corresponding operators for the particle confined in 
the one-dimensional infinite square well [9], and are also connected with the two-dimensional 
harmonic oscillator in a sector operators [16]. 

The reader may also examine the actions of the corresponding operators in the j direction 
recognizing that they lead to the companion cosine functions satisfying the Neumann boundary 
condition instead of the Dirichlet one of equations (8), (9). 

For a fixed value of p, the radial dependent operators in the combinations 

d p 
dp p 

guarantee the removal of the singularities at p — > and p — > oo, including the annihilation of the 
eigenstates p^e~ p I 2 when the upper signs are chosen. For the states with radial excitations, we 
adapt the recurrence relations for the confluent hypergeometric functions (13.4.10) and (13.4.11) 
in [1]: 



o + z~ ) M(a,b,z) = aM(a + l,b,z), 
dz ' 



b-a- z + z— ) M(a, b, z) = (b - a)M(a - 1, b, z) 



(25) 
(26) 



with the correspondences a = — n p , b = p + 1 and z = p 2 in equation (13). The first one with 
a + 1 = —(rip — l)is the lowering operator, while the second one with a — 1 = — (n p + 1) is the 
raising operator. The interested readers can check these actions between the consecutive entries 
in Table 1. 

Since the harmonic oscillator in the z-direction is free, the actions of the corresponding 
operators in equation (19) on the Hermite polynomials of Table 2 are the familiar ones. 

In the spherical coordinates, the raising and lowering operators in the ip angle are those in 
equations (23), (24), while those in the radial coordinate are also given by equations (25), (26) 
with the correspondences a = —n r , b = A + |, and z = r 2 in equation (16). Their actions can be 
verified by applying them in the entries of Table 3. For the polar angle the recurrence relations 
for the associated Legendre functions (8.5.3), (8.5.4) in [1] 



Xz 



1) 



d_ 
dz 



P^z) = (X + p)P^ 



(A + l)z+(A 2 -l) 



dz 



P^z) = {\-p + l)P^ +1 (z) 



(27) 
(28) 



can not be applied directly to the eigenfunctions in equation (17) due to some of the remarks 
in that paragraph. 



Ladder Operators for Quantum Systems Confined by Dihedral Angles 



7 



Nevertheless, equations (27), (28) suggest the structure of the correct lowering and raising 
operators for the polar angle eigenfunctions with parity ir = (—) ne - 



d 

dz 



2 Fi -no, 2// + n + 1; fi + 1; 



1-z 



7r=(-l) n 9 



2 Fi -(no - 1), 2/i + n e ; [i + 1; 



1 - z 



7r=(-l)"fl 



(29) 



az 



2 F X -(n e + l),2fi + n e + 2; n + 



2 Fi -n e , 2/i + n e + 1; /i, + 1; 



1 - z 



*=(-)•* 



1 - z 



7r= (_)"8+l 



(30) 



Their applications to the entries in Table 4 lead by mathematical induction to the general 
results: 



jy( n e) _ -y £>{ne) 



1 

n 9 ' 



a 



(n e ) 



1 

no' 



A 



(ne) 



1, B 



and in particular, = for annihilation of the eigenfunction with no = 0, which is 1 and also 



( n e) 



1 



2/i + n e + 1 



1 



2/i + n e + 1 



(31) 
(32) 



1, and B 



(0) 



has a vanishing derivative, making the values of B^ and ci°^ irrelevant; also A ? 

and irrelevant. Notice also the equal magnitudes and opposite signs of the respective B 
and C coefficients leading to the common factors =f(z 2 — 1) of the derivative terms in equa- 
tions (29) and (30), respectively. The values of A^ B ^ and B^ 9 ^ guarantee the vanishing of the 

term with the no power of z, while the values of A^ e ^ and B^ B ^ guarantee the highest power 
(no + 1) of z and its correct coefficient in the r.h.s. eigenfunction in equation (30). 

The energy eigenvalues in equations (15) and (18) contain the same zero-point, radial and 
axial or polar angle contributions as in the free harmonic oscillator. The difference for the 
oscillator confined by dihedral angles corresponds to the \x contribution which according to 
equation (12) is quantized with integer values of and spacings inversely proportional to the 
confining angle (po. 

Let us start with the case of tpo = it, for which the dihedral angle becomes a plane and fi 
becomes n«; consequently, its contributions to the energy are on the same footing as those of the 
other degrees of freedom. Next, let us consider the confining angles in the interval it ^ tpo ^ 2ir 
for which the /i energy contribution is less than those for the other degrees of freedom, taking its 
lowest value of one half for ipo = 2tt; this is the quasi-free limit for the harmonic oscillator, being 
excluded from the xz half plane (0 < x < oo, y = 0, — oo < z < oo), and definitely different from 
the free harmonic oscillator. On the other hand, for the interval tt ^ ipo > 0, as the dihedral 
angles is closing starting from the plane, the /i contribution to the energy becomes larger than 
the contributions from the other degrees of freedom, tending to infinity as ipo — > when the 
moving meridian half-plane of the confining angle approaches the position of the xz half plane. 



4 Interbasis expansions of degenerate eigenfunctions 
in circular cylindrical and spherical coordinates 

In this section the complete eigenfunctions of the harmonic oscillator confined by dihedral angles 
in circular cylindrical and spherical coordinates are compared, illustrating explicitly the rela- 
tionships of their sets with the same energy, as counterparts of [8, 12, 14] for the free oscillator. 
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The complete eigenfunctions are the products of the common factor of equation (12), and 
the respective companions of equations (13), (14) and equations (16), (17): 

^n pt m z (fi, V, z) = R npP (p)$ p (ip)Z nz (z), ip nrngtl (r,6,ip) = flnmj^r)^^^^). 

For a chosen value of p, the eigenstates without excitations in the other degrees of freedom 
are the same in both coordinate systems, which is expressed as 

\n p = 0, p, n z = 0) = \n r = 0,ng = 0, p) 

using the ket notation. The eigenstates with excitations in the other degrees of freedom share 
the common factors in equations (13), (14) and (16), (17) as already discussed in the paragraph 
after equation (18). Correspondingly, their comparison can be restricted to that of the products 
of the polynomials in Tables 1 and 2 and in Tables 3 and 4, respectively. Such a comparison will 
be implemented for the eigenstates involving the successive rows in those tables, and making 
use of the connections between the respective coordinates, equation (7). 

The lowest excitation eigenstates, with negative z-parity and energy (// + 5/2) are also the 
same: 

\n p = 0, p, n z = 1) = \n r = 0, n g = 1, p) 
since z = r cos 9. 

The next sets of excited states with positive z-parity and energy (fi + ^), \n p = l,p,n z = 0), 
\n p = 0, p, n z = 2) and \n r = 1, ng = 0, p), \n r = 0,ng = 2, fi) involve the products of the respec- 
tive polynomials: 



1 2 \ 2 (2m + 3) cos 2 9-1 

i -r • 1, r 



V M+l J ' 2 (/ u + 1 ) 

By using r 2 cos 2 6 = z 2 and r 2 = p 2 +z 2 , the last pair can be written as the linear combinations 
of the first pair: 

1 _W + ^\ = V£±»( l _J?_\ 1 _(2 2 2 _i) (33 ) 

2^ + 3/ 2/x + 3 V /i + l/2^ + 3 v ; ' vy 

( 2 , + 3)^-^-^)_ 1 ^^Vfcl) (!j! _, 



2(/i + 1) 2p + 3 V A* + 1 / 2/i + 3 

The relationships between the normalized eigenfunctions follows from these relationships 
between the polynomials 

K = i, n = o, p)\ = _ y-r:\ /|n p = i,p,n z = oy 

\n r = 0, n g = 2, //)/ \J 2^f3 + \/ 2 2^+l/ \\n p = 0, /i, n 2 



In the same way we have also obtained the following results for the successive eigenstates 

9 m + I! „+ 13 



with higher excitations and energies E = p + ^,p + ^,p + ^§, respectively, 



\n r = l,n = l,p)\ /^g) /K = l,p,n z = l) 

\n r = 0,n$ = 3,p)J \J 27^5 V 2^+5 J \\n p = 0, p,n z -- 
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/|20m>\ 
I 12m) 

VloW 

/|21 At )\ 

1 13 M ) 
V|05/i)/ 



/ / 4( At +l)(M+2)" 
(2^+3) (2^+5) 

4(/x+2) 
(2 M +3)(2^+7) 



(2 M +5)(2/i+7) 

4(^+l)( M +2) 
(2^+5) (2^+7) 
12(^+2) ~ 
(2/H-5) (2/4+9) 

15 



4(^+1) 
(2/i+3)(2^+5) 

2^+1 
V (2 M +3)(2 M +7) 
12(/j+2) 
(2/H-5)(2/i+7) 

12(^+1) 
(2/x+5)(2/H-7) 
2/t-l 
V (2 M +5)(2 M +9) 
20(^+2) 



(2^+3) (2/i+5) 
12(^+1) 
(2/i+3)(2/x+7) 

4( M +l)(^+2)~ 
(2 M +5)(2^+7) / 



15 



\Y (2 M +7)(2 M +9) V (2 M +7)(2 At +9) 



(2/H-5)(2/j+7) 
20(^+1) 
(2^+5) (2^+9) 

4( M +l)(M+2) 
(2 M +7)(2^+9) / 



/|2 M 0)\ 

|l/i2) 

V|0/x4)7 
/|2/il)\ 

|l/x3) 

V|o M 5>y 



Notice the orthonormality of the transformation matrices, implying the Solomonic distri- 
bution of the coefficients in equations (33), (34) and their counterparts between the entries 
in the matrices and the normalization constants of the respective eigenfunctions. The sets 
of eigenstates come by pairs of singlets, doublets, triplets, and in general multiplets of order 
(2n p + n% + 2)/2 = (2n p + n z + l)/2 with n° z = n e z + 1, etc., of positive and negative z-parity, 
and with equal successive energy spacings of Huj, for fixed values of \x. The relationships in this 
section are valid for any chosen values of ii from equation (12), for any excitation and for 
each angle of confinement <p$ . 



5 Eigenstates and eigenenergies of the hydrogen atom 
confined by dihedral angles 

The Schrodinger equation for the hydrogen atom is separable in spherical coordinates, equa- 
tions (7), as well as in parabolic and prolate spheroidal coordinates [2, 3, 4, 6, 10, 13, 21, 22, 24]: 



x = ^rj cos p = / ' \J (u 2 — 1)(1 — v 2 ) cos ip, y = £77 sin p = fy/ (u 2 — 1)(1 — v 2 ) sin ip, 

fuv. (35) 



e-r, 2 



The eigenfunctions of the hydrogen atom confined by dihedral angles in the respective coor- 
dinates must satisfy the equations (9) and their counterparts 



^(£,77,^ = 0) = 0, v p (^,v? = v?o) = o, 



*p PS (u, v,<p = 0) = 0, V(u, v,cp = ip ) = 0. 



PS 1 



(36) 
(37) 



They also share the eigenvalue equation and its solutions in equations (10), (12). In spherical 
coordinates, the eigenvalue equation for the square of the angular momentum is also the same 
as in Section 2, and for the spherical radial Hamiltonian: 



2m P 



1 d 2 d A(A + 1) 
r 2 dr dr r 2 



e 
r 



R{r) 



ER (r) 



In parabolic coordinates, the z-component of the Runge-Lenz vector is the constant of motion 
leading to the separated ordinary differential equations: 



2m e 

h 2 



1 d d fjr 
Zdt dt ~ Z 2 
1 d d 



V- 



/i- 



2m e \ 77 dr] dr] 7] 



Erf 



E (0 = A^E (0 , 
H(r ] ) = A v H(r l ), 



(38) 
(39) 
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where the separation constants must satisfy the condition 

A). + A v = 2e 2 . (40) 

For the bound states with negative energy of the hydrogen atom, these equations are equivalent 
to that of the two-dimensional harmonic oscillator in two dimensions in the circular radial 
coordinate, equation (11). 

The constant of motion in prolate spheroidal coordinates is a linear combination of the square 
of the angular momentum, the z-component of the Runge-Lenz vector and the Hamiltonian [15], 
leading to the separated ordinary differential equations 



du du u 2 — 1 



2me 2 fu 2me z Ef 2 u 2 



h 2 h 2 



j I" d , 2 \ d V 2 1 2me 2 fv 2me 2 Ef 2 v 2 1 
\[dv^ )dv~T^j 2 \ + H 2 + ~ ~H 2 J 



U (u) = -AU (u) 
V(v) = AV(v) . 



These two equations are the same in the respective domains 1 ^ u < oo ^ and — 1 ^ v ^ 1. 

Their remaining radial and angular solutions have the same forms as those of the free hydrogen 
with the substitution of the integer magnetic quantum number m with the eigenvalue fi in 
equation (12). 

In spherical coordinates, the polar angle eigenfunction of the square angular momentum is 
the one in equation (17), and the radial eigenfunction of the energy becomes 

<.a(0 = N nrX r x e- r /» a »M(-n r , 2A + 2, 2r/va ), (41) 

with A = n$ + (j,. The energy eigenvalue E = —e 2 /2aQV 2 involves the non-integer principal 
quantum label 

v = n r + rig + n + 1. (42) 

In parabolic coordinates, by adopting the harmonic oscillator parametrization in equa- 
tions (11) and (15) 

1 e 2 

E = - -moj 2 = - — , A% = fvjj(2n^ + /i + 1), A v = Huj(2n ri + /i + 1), 

the identification of their frequency ui = e 2 /hv, length unit H/moj = uao, and the connection 
between the quantum labels follow immediately: 

v = ri£ + n v + n + 1, (43) 

as the counterpart of equation (42), and consistent with equation (40). 
The respective eigenfunctions share the same forms 

= N n ^e^ e/2uao M(-n^fi + 1, f/va ) , (44) 
H n ^ = N nri ^e-^/ 2 " ao M(-n v ,^ + 1, V 2 /va ) , (45) 

which coincide with that of equation (13) for the harmonic oscillator in the circular radial 
coordinate. Notice also that £r/ in equations (35) coincides with p = rsva.9 in equations (7), 
and that (£ 2 + rj 2 )/2 = r using equations (35) ensuring that the products of the factors in 
equations (44), (45) and equations (41), (17) share the same factors removing the singularities 
at the z-axis and at infinity. 
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Table 5. Confluent hypergeometric polynomials. 



M(-n r ,2X + 2,2r) 



1 

2(2r) 



2r 
2A+2 



3(2r) 



+ 



2A+2 
3(2r) 2 



+ 



(2r) 2 
(2A+2)(2A+3) 



(2r) 3 



2A+2 (2A+2)(2A+3) (2A+2)(2A+3)(2A+4) 



In prolate spheroidal coordinates the spheroidal and hyperboloidal eigenfunctions also share 
the same differential equation and forms: 



U nu ,(u) = N null (u 2 - l)^V/"/-°S£>), 
W) = ^(1 - ^V^S^), 
with polynomials 



on 



E c «(^- 1 ) 



(46) 
(47) 



(46 



s=0 



s=0 



with common coefficients satisfying three-term recurrence relations [6, 15] 



2/ 



(n u + n„ + 1 - s)c s _i + 



s ^ s + 2/i + 1 
+ 2(s + l)( s + /i + l) Cs+ i = 0, 



2/ 
a ^ 



.4 



(49) 



and in different ranges of their variables l^u<oo, — l^v^l. Now, equations (35) 
and (7) lead to the identification of f \f (u 2 — 1)(1 — v 2 ) = n sin^i and f(u + v) = r±, where the 
subindex 1 indicates that the spherical coordinates are referred to the position of the nucleus 
as the origin and located at the focus (x = 0, y = 0, z = — /). Correspondingly, the pro- 
ducts of equations (46), (47) and equations (41), (17) also share the same factors removing the 
singularities at the z-axis and at infinity. 
The principal quantum label becomes 



n u + n v + fi + 1. 



(50) 



Comparison of equations (42), (43) and (50) with a common value of [i lead us to recognize 
the common value of their energy for the states with 

n r + n e = ri£ + n v = n u + n v , 

and with a degeneracy one unit above such sums. 

Table 5 illustrates the explicit expressions of the polynomial confluent hypergeometric func- 
tions of equation (41). Their polar angle companions are those in Table 4. 

The parabolic coordinate confluent hypergeometric polynomials of equations (44), (45) coin- 
cide with those in Table 1 with the respective replacement of p with £ and rj, as already discussed. 

Table 6 illustrates directly the ratio of normalization constants and polynomial factors for 
excited states in spherical and parabolic coordinates, respectively. 

The products of the spheroidal polynomials in equations (48) for the excited states with 
n i + n 2 = 1)2, obtained from the recurrence relation of equation (49), have the same forms 



1 + 



Aj(u-iy 

2(1 + /x) . 



1 + 



Mv - 1) 
2(1 + //) 



(51) 
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Table 6. The ratios of normalization constants iV ni „ 2(tt / /Nqq^ and polynomial factors for excited states. 



States 



Spherical coordinates 



Parabolic coordinates 



01/x 
10// 
02 fi 

20/i 



r cos 8 



/2+fi a u 
2(1+M) f 1 



/2+Ji 



4(1+M) 



2 [-l+(3+2/i) cos 2 6] 



(2+ M )(3+/i)(3+2/i) 
\f$,r cos 8 ( 



(a v) z (l+At) 



aov(2+fi) 



3+fJ. 



X 1 



2r 



2r 2 



a v(l+p) agi/ 2 (l+M)(3+2^) 




(1+M)y/2(1+^) 



for the states \n u = 0,n v = 1, /x) and |n u = 1, n v = 0, /i) where with 2 = 1,2 are the ordered 
roots of the quadratic equation 



A 



2 i,.,,, A L\ A _^L 



2/i + 2 



a fi 



0. 



(52) 



Similarly, the states \n u = 0,n v = 2,/i), \n u = l,n v = 1,/x) and \n u = 2,n v = 0, /i) share the 
products 



1 + 



Mu-V) _ [A- 2(1 + M )](4/ + - 1} 

2(1 + 



1 + 



8(l + fi)(2 + n)v 
[Ai - 2(1 + fi)](4f + Aiv)(v - l) 2 



2(1 + /i) 8(l + /x)(2 + /i)i/ 

where Ai with i = 1,2,3 are the ordered roots of the cubic equation 



(53) 



A 3 + (-8-6/J+ — )A 2 + 



+ 



a v J 

16/(1 + ^(2/1 + 3) 64/ 2 (/i + l) 



12 + 20. + 8^- 16/ ^ + 4) + ^ 



(a /i) 2 



A 



a /i 



(a /i) ; 



0. 



(54) 



The generalization for any fixed value of n u +n v involves the tridiagonal matrix representation 
of equation (49) of dimension (n u +n v + l) x (n u +n v + l). The corresponding secular determinant 
generalizing equations (52) and (54) is an algebraic equation for the eigenvalues Ai of degree 
n u + n v + 1. The diagonalization of the matrix yields the numerical values of A^s and the 
coefficients c s (Ai) for the respective polynomials in equations (48). 



6 Raising and lowering operators 

On the basis of the equivalence of the hydrogen atom eigenfunctions in parabolic coordinates, 
equations (44), (45), and the harmonic oscillator eigenfunctions in the circular radial coordinates, 
equation (13), they also share the same kind of raising and lowering operators of equations (25) 
and (26) with the correspondences a = — ng, — n rv b = /x + 1, z = £ 2 /vao£, rj 2 /vclq, 

^/i + 1 + n 5 - z + z-^j M (-ri£, n + 1, z) = (p + 1 + n € ) M (- (n$ + 1) , p + 1, z) , 
(~ n? + z Tz) M ^~ n ^ ^ + X ' z ) = ~ n « M (- ( n H ~ 1) > M + 1, ^) . 
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In spherical coordinates, both the hydrogen atom and the harmonic oscillator share the same 
polar angle eigenfunctions of the square of the angular momentum and z-parity, equation (17), as 
well as their raising and lowering operators, in equations (29)-(32). For the hydrogen atom radial 
eigenfunctions of equation (41), the corresponding operators also follow from equations (25), (26) 
using the correspondences a = —n r , b = 2A + 2, z = 2r jva^: 



2A + 2 + n r - z + z—\ M(-n r , 2A + 2, z) = (2A + 2 + n r )M(-(n r + 1)) + 2(2A + 2, z), 



The interested readers can easily check their connecting actions between the successive entries 
in Table 5, including their extension. 

Concerning the raising and lowering operators for the spheroidal polynomials of equa- 
tions (48), the inspection of their explicit forms for the lowest excited states in equations (51)- 
(54) lead us to their respective identifications. In fact, for each fixed value of n u + n v , the 
n u + n v + 1 degenerate eigenfunctions share the same form differing in the associated eigen- 
value Aj. Consequently, the successive changes in the neighboring eigenvalues A\ — > A2 — > 
A 3 -»■ > A nu+nv -> A nu +n v +\ translates into the raisings for the spheroidal quantum num- 
ber n u = n u = 1 — > n u = 2 — y ■ ■ ■ — > n u = n Umax - 1 — > n u = n Umax with the corresponding 
lowerings for the hyperboloidal quantum number n„ max — > n 1)max — 1 — > n Viaxx — 2 n« = 

1 ->• n v = 0, in such way that n Umax = n„ max = n u + n v . 

On the other hand, for the polynomial solutions in equations (48), their common expan- 
sion coefficients depend on the respective eigenvalue c s (A nu ) where we replace the index i = 
1,2,..., iV max + 1 ordering the eigenvalues with the quantum number n„ = 0, 1, 2, ... , iV max 
counting the number of spheroidal nodes. Notice that both polynomials have the same form 
and are of degree n u + n v , their difference being associated with the distribution of their nodes 
in the spheroidal and hyperboloidal degrees of freedom, in a complementary way. Their total 
number is n u + n v + 1 associated with states of the same energy. The powers of u — 1 or v — 1 
in which they are expanded appear in the same number in equations (48). Both finite sets of 
linearly independent sets of power functions can be expressed as linear superpositions of the cor- 
responding sets S nu (u) and S nv (v), via the inverse transformation based on the inverse matrix 
of the coefficient c s (A nu ). These relationships have their counterparts for the familiar special 
functions, in polynomial and infinite series forms. 

7 Interbasis expansions of degenerate eigenfunctions 

in spherical, parabolic and prolate spheroidal coordinates 

Following the same order as in Section 4, the normalized eigenstates without excitations in the 
non-common degrees of freedom are the same in the three coordinates systems: 

\n r = 0, uq = 0, ijl) s = |ri£ = 0, n v = 0, ji) P = \n u = 0, n v = 0, n) PS ■ 

From the entries in Table 6, the connections between the normalized eigenstates in spherical 
and parabolic coordinates become: 




n r M(-(n r - 1),2A + 2, z). 
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and 



/|02/x) s \ 
111/") 5 



i+jj, 

2(3+2^) 
1 
V2 



3+2^ 





/ 2+M / 1+M / 

\Y 2(3+2 M ) V 3+2 M V 



2^ (l°^ P \ 
|1^) P 

\\20f,) p J 



l 

V2 



2+A* , 
2(3+2^)/ 



where the orthonormality of the transformation matrices is obvious. These are the counterparts 
of the interbasis expansions for the free hydrogen atom [20]. 

Similarly, from the spherical entries in Table (6) and equations (48) for the prolate spheroidal 
polynomials, the transformation matrices between the respective products of polynomials are 



2(1+M) 
f(A 2 -A 1 ) 



4L 



2(1+/*) 



2(1+M) 
7(^2 -Aa) 



2(1+^)(A 2 -Ai) 



A, 



A 2 



2(l+/i)(A 2 -Ai) 



and 



8(l+/ t ){2(l+M)[4+3M-(A 2 +A 3 )]+A 2J 4 3 } 
(Ai-A 2 )(Ai-A 3 ) 
8/(l+At)[A 2 -2(l+ A t)][A 3 -2(l+ M )] 
(Ai - A 2 ) (Ai - A 3 ) ( Ai+2/) [Ai -2(1+/Lt)] 

16/ 2 (2+ M )[A 2 A 3 +4/(1+m)] 



8(l+/i){2(l+/ J )[4+3M-(A 3 +A 1 )]+A 3 Ai} 
(Ai-A 2 )(Ai-A 3 ) 
8/(l+^)[A 3 -2(l+^)][Ai-2(l+ M )] 
(Ai-A 2 )(Ai-A 3 )(Ai+2/)[Ai-2(l+/i)] 
16/ 2 (2+ m )[A 1 A 3 +4/(1+m)] 



( Ai - A 2 ) ( Ai - A 3 ) ( Ai +2/) (3+2/1) [Ai -2(1+ M )] (A 1 -A 2 )(A 2 -A 3 )(A 2 +2/)(3+2 M )[A 2 -2(l+ At )] 

8(1+m){2(1+ m )[4+3m-(A 3 +Ai)]+A 3 Ai} 
(Ai-A 2 )(Ai-A 3 ) 
8/(l+/x)[Ai-2(l+/x)][A 2 -2(l+At)] 
(Ai-A 2 )(A 3 -A 2 )(A 3 +2/)[A 3 -2(l+ M )] 

16/ 2 (2+ M )[AiA 2 +4/(l+ At )] 



(A 1 -A 3 )(A 2 -A 3 )(A 3 +2/)(3+2 M )[A 3 -2(l+ M )]_ 

These interbasis expansions are the counterparts of those for the free hydrogen atom [17, 18]. 



8 Discussion 

The difference between the free harmonic oscillator and hydrogen atom and their confined by 
dihedral angle counterparts resides in the boundary conditions of equations (8), (9) and (36), (37) 
in the respective coordinate systems. The latter lead to the eigenfunctions and eigenvalues of 
the square of the z-component of angular momentum of equation (12), where \i replaces the 
familiar integer magnetic quantum number m. 

Nevertheless, both the free and the confined- by-dihedral-angle quantum systems share their 
superintegrability and exact solvability. Admittedly, the separability and integrability is re- 
stricted to the coordinate systems sharing the dihedral angle coordinate: circular cylindrical 
and spherical for the harmonic oscillator, and spherical, parabolic and prolate spheroidal for the 
hydrogen atom. Sections 3 and 6 illustrated the identification and construction of the raising 
and lowering operators for the eigenfunctions of l? z , the transverse and axial contributions to the 
Hamiltonian of the isotropic harmonic oscillator, its complete spherical radial Hamiltonian, the 
square of the angular momentum and z-parity operators; and their extensions and adaptations 
for the Hydrogen atom in the respective coordinates. These operators have their free harmonic 
oscillator and Coulomb potential counterparts in [5]. Concerning special functions, the reader 
is invited to recognize the roles of the fractional Fourier basis, equation (12), and the square 
of the angular momentum and z-parity eigenfunctions equation (17). Additionally, Sections 4 
and 7 illustrate the relationships between the complete degenerate eigenfunctions of the confined 
quantum systems in pairs of the alternative coordinate systems. 
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Some of the eigenfunctions and their raising and lowering operators may also be adapted 
to other quantum or electromagnetic systems sharing the same constants of motion and con- 
finement situations. As specific examples the free particle confined by dihedral angles, and 
nondiffracting vortex beams [11], can be mentioned. 
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